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Abstract 

We investigate several fluctuation effects in high-energy hadronic and nuclear 
collisions through the analysis of different observables. To introduce fluctu- 
ations in the initial stage of collisions, we use the Interacting Gluon Model 
(IGM) modified by the inclusion of the impact parameter. The inelasticity 
and leading-particle distributions follow directly from this model. The fluc- 
tuation effects on rapidity distributions are then studied by using Landau's 
Hydrodynamic Model in one dimension. To investigate further the effects 
of the multiplicity fluctuation, we use the Longitudinal Phase-Space Model, 
with the multiplicity distribution calculated within the hydrodynamic model, 
and the initial conditions given by the IGM. Forward-backward correlation is 
obtained in this way. 
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I. INTRODUCTION 



One of the main characteristics of the high-energy hadronic or nuclear collisions is the 
existence of large event-by-event fluctuations: fluctuation of multiplicity, of particle species, 
of inelasticity, of momentum distribution, of impact parameter, and so on. Usually, one tries 
to describe an average trend of such a phenomenon, of multiparticle production, by using 
for example hydrodynamic models HQ and with success. However, if one analyzes the data 
more closely, one observes that in a given experimental setup, even under the same initial 
conditions of colliding objects, events with different final state configurations take place. This 
fluctuation has either a quantum mechanical or a statistical origin or even simply associated 
with the impact parameter. The so-called inclusive data for the final-particle distributions 
are the averages over such event-by-event fluctuations for a given set of experimental initial 
conditions. In the usual application of hydrodynamic models, to describing the inclusive 
data we presumably expect, by means of a sort of ergodic assumption ||, that the average 
over event-by-event fluctuations is replaced by a statistical ensemble average. However, not 
all the average of physical fluctuations can be expressed in terms of the above average over 
statistical ensemble of the constituent configurations. For example, the impact-parameter 
and quantum-mechanical fluctuations that occur in the initial condition of each event can 
never be averaged out with the use of the ergodic hypothesis. Besides, a description in 
terms of the average quantities is clearly not satisfactory in treating such quantities as 
multiplicity, inelasticity, semi-inclusive rapidity distributions, correlations among particles, 
etc. The main aim of this work is to discuss the effects of such fluctuations on the observed 
quantities and to present our attempt to include them in the description of the data. 

In the following, we shall give in the next two sections a brief account of a modified version 
011 of the Interacting Gluon Model (IGM) |J, which will be used to generate the event- 



1 The field-theoretical foundation of hydrodynamic model has first been given by M. Namiki and 
C. Iso § 



dependent fireballs. In Section II, we discuss how the impact-parameter fluctuation can be 
treated and in Section III, how IGM incorporates the energy and momentum fluctuations 
of the partons to obtain the event- dependent energy and momentum of the central fireball. 
Comparisons with some data, showing the effects of these are also performed there. These 
fluctuations are nevertheless not enough, for there are many data which require consideration 
of other kinds of fluctuations. In Section IV, we describe how the multiplicity fluctuation 
may be treated, concentrating especially in pp collisions and their multiplicity distributions. 
To treating the energy and momentum fluctuation of the observed particles, once the mass 
of the fireball and its multiplicity are defined, we use the longitudinal phase-space model. 
This is discussed in Section V. It is shown that these two kinds of fluctuations are essential 
in describing the forward-backward correlations. Conclusions are summarized in Section VI. 



II. IMPACT-PARAMETER FLUCTUATION 

In any collision process between particles or nuclei, the impact parameter cannot be fixed 
a priori. This is in part due to the quantum mechanical uncertainty, but, even if we could 
theoretically define the trajectories of the incident particles like in heavy-ion collisions where 
the incident objects are nearly classical, it would equally not be possible in practice due to the 
actual experimental conditions. We may recall that there exist some experimental techniques 
to discriminate, a posteriori, the central from the peripheral collisions in such reactions. But 
they do not eliminate ponderable uncertainties. So in any realistic description of hadronic 
or nuclear collisions the impact-parameter fluctuation must be included. In previous works 
we have studied this fluctuation in connection with IGM and hydrodynamic model 
and shown that it affects the observables such as the inelasticity, leading-particle spectra, 
rapidity distributions of produced particles in a significant amount. 

The impact parameter b defines, in the first place, the probability density of occurrence 
of a reaction (apart from the normalization) F(b) = 1 — |5(6)| 2 , where the eikonal function 
is written as 



\S(b)\ 2 = exp{-cJdb' J db ,/ D A 0)D B (b'')f(b + b'-b'')} = exp{-Ch AB (b)}, (1) 

where D A ^{b) are the thickness functions of the incident particle (nucleus) A and B re- 
spectively. C is an energy- dependent parameter and, though not necessarily so, in studying 
pp, p-nucleus and nucleus-nucleus collisions, we assume it to be universal, so that it may be 
determined by the condition / F pp (b) db = a % pp l (\/s) for pp collision. Notice that, because 
of this, once pp cross-section is fixed, the AB cross-section cr^J (\/s) = J F A s{b) db may be 
calculated by using (|J). We take, as an input 0, 

a l p n p el = 56 (v^r L12 + 18.16 (v^) ' 16 - (2) 

— # — * — # 

The function f(b) in ([]]), subject to the constraint / f(b)db = 1, accounts for the finite 
effective parton interaction range (with the screening effect taken into account). The simplest 
choice of fib) would be the point interaction 5(b), but we prefer to parametrize it as a 
Gaussian with a range ~ 0.8 fm, which is more consistent with the character of the strong 
interaction and also describes better the data. For proton, we parametrize D p (b) as a 
Gaussian distribution. So, we have eventually D p (b) = f(b) = (a/V) exp(— ab 2 ), with 
a = 3/(2 R 2 ) , where Rp ~ 0.8 fm is the proton radius. For nuclei, we take 

D A {b)= [ p A {b,z)dz= f — — dz , (3) 

J-oo J-oo 1 + exp [(r — R ) /d\ 

where Rq = tqA 1 ^ 3 , r = 1.2 fm, d = 0.54 fm and p^(r) is normalized to A. Thus, in the 
particular case of pA collisions, we get 

h pA (b) = a [°° db'b'D A (b') Io(abb') e - a{b2+b ' 2)/2 , (4) 

where Iq is a modified Bessel function. 

Besides, the impact parameter determines the size of the fireball, because as b increases 
the overlap of the hadronic matter becomes smaller and consequently so does the average 
mass of the fireball generated in the collision. We incorporate this effect by writing the 
parton momentum distribution functions as 



G A (x,b)=D A {b)/x, G B (y,b)=D B (b)/y, (5) 

where x and y are the Feynman variables of partons in A and B, respectively, in the equal- 
velocity (e.v.) frame. With this, we are assuming that the parton momentum distribution is 
independent of the particular type of nucleus and also equal for whole the nucleus, the only 
difference being their density (where the thickness D(b) is large, more partons of a given 
momentum are found). Here, it is presumed that the same physics describes pp, pA and 
AB collisions and that no correlation exists among the nucleons inside each of the incident 
nuclei. Then, given an impact parameter 6, the density of parton pairs with momenta 
(x, —y)\fs/2 that fuse contributing to the fireball formation is written as 

w(x, y-b) = J db' J db" G A (x, b") G B (y, b") a gg (x, y) f(b+ b' - b") 9 (xy - M^ in / s) 

= h AB (b)w(x,y) , (6) 

with 

w(x,y) = [o- gg (x,y) /xy] 6 (xy - M^ in /s) , (7) 

where M min = 2m n and the parton-parton cross-section is parametrized asQ a gg (x,y) = 
a I (xys). Observe that in (||), b dependence is factorized out. 

III. ENERGY AND MOMENTUM FLUCTUATION OF THE CENTRAL 
FIREBALL: INTERACTING-GLUON MODEL 

Interacting Gluon Model (IGM) |J is a simple QCD motivated model, especially de- 
signed to create the initial conditions for hydrodynamic descriptions, by incorporating in 
an intuitive way the microscopic fluctuations in the initial stage of the collision. It is based 
on an idea that in high-energy collisions valence quarks weakly interact so that they 

2 In following the original IGM ||, we had parametrized (0) with a™ el in the denominator. In 
|J, we have redefined a as an adimensional constant like in this expression. 



almost pass through, whereas gluons interact strongly, producing an indefinite number of 
mini-fireballs, which eventually form a unique large central fireball (all possible qq sea quarks 
are, in this model, "converted" to equivalent gluons). One of the nice features of this model 
is its easy handling. However, the main drawback of the original version was the neglect 
of the impact-parameter. In [[|,[5|, we have improved it, by including the impact-parameter 
fluctuation as explained in the previous section. By reinterpreting the partons in that sec- 
tion as gluons, since only these are assumed to interact, the probability density x(E, P; b) 
of forming a fireball with energy E and momentum P at a fixed b is obtained as follows. 
We assume that the colliding objects form a fireball, via gluon exchanges, depositing in it 
momenta x(b)y/s/2 and —y(b)y/s/2, respectively. Let be the number of gluon pairs that 
carry momenta Xiy/s/2 and —yiyfs/2. Thus, 



E 



TLjXi 



x(b) and = y(b) . 



(8) 



In what follows, we will omit the explicit b dependence of x and y in order not to overload 
the notation. The energy and momentum of the central fireball in the e.v. frame of the 
incident particles are given by 



E = (x + y) 



[X 



y)\Ts/2 



(9) 



and its invariant mass M and rapidity Y are respectively 



M = ^Jsxy = ka/s and Y = (1/2) In (x/y) . 



(10) 



With these notations, we can follow the prescription given in || and write the relative 
probability of forming a fireball with a specific energy and momentum as 



r(x, y\ b) ~ exp{-X J G X}/ 



7T 



detfG^ 



where 



X 



x—{x) 



y-(y) 



G 



(x 2 ) (xy) ^ 
(xy) (y 2 ) 



(11) 



with the notation 

( x m y n ) = Jdx' J dy'x' m y' n w(x',y';b), (12) 

or, in terms of E and P, 

T(E, P; b) ~ exp{-ai(E - {E)f - a 2 P 2 }, (13) 

where oi = [s((x 2 ) + (xy))]" 1 , a 2 = [s((x 2 ) - {xy})}' 1 and (£) = ((x) + (y))y/s/2 
(don't confuse this notation with the average value; it is not because w(x,y;b) is not nor- 



malized). Apparently, T(E, P;b) in (13) is normalized. However, both E and P are bounded 
because of the energy- momentum conservation constraint. It is also constrained by M > 
M m i n = 2m n . So, we put some additional factor Xo(b), 

X (E,P;b) = Xo(b)T(E,P;b), (14) 

such that 

Up (dE X (E, P; b) 6(VW=P* - M min ) = . (15) 

J J (Tad 



As implied by @ and remarked there, the gluon momentum distribution is independent 
of the particular type of nucleus, the only difference being their density. So, in the integral 



(12]), x' and y' vary from some lower limit, defined by a/ sx'y' = M min , up to 1, corresponding 
to the complete neglect of any collective effect of the nucleons in a nucleus. On the other 
hand, the integration limits of (|15|) are chosen differently, x and y in ([8]) may be larger than 



1, because gluons from different nucleons may contribute to give the fireball a momentum 
transfer that is larger than \fsj2 , which is just the incident momentum of a single nucleon 
in our e.v. frame. We take as the upper limit of x and y the overlap h,A,B{b), whenever it is 
larger than 1, and limited to A and B, respectively (x or y remains < 1 , if it corresponds 
to a proton). When hA,B{b) < 1, we take it = 1, because in such a case just a single nucleon 
of each nucleus interacts. 



A. Inelasticity distributions 



The concept of inelasticity, understood as the fraction of the incident energy E which is 
lost while a particle interact with another one, is crucial in cosmic-ray data analysis where the 
primary mass composition, hence an important information about the Universe, is deduced 
by using models of cascade development with appropriate inelasticity distributions^ and 
cross sections as the inputs. In this case, the usual definition is k = (Eq — E') / Eq , where 
E' is the leading (or surviving) particle energy. At high energy, a nucleus hardly supports 
any collision without suffering a breakup, so the term inelasticity in this sense is commonly 
reserved to a hadron projectile and, in this paper, we shall restrict the consideration only to 
proton incident on nuclear targets. However, in a wider sense, as the fraction of the incident 
energy Eq which is materialized into produced particles, it is also important in connection 



with the production of a quark-gluon plasma in heavy- ion collisions [p3|HT4|6]1 ■ In Ref. ||, 
n appearing in ( |T0| ) is called inelasticity and this is the quantity of interest in heavy-ion 
collisions. In the present paper, we shall adopt the usual definition k given above and, as 
for k, call it simply k. In any case, it is clear that we can talk about inelasticity distribution 
only when there exists some event-by-event fluctuation. 

Having obtained x(E, P; b) , we can readily compute such distributions. The k- 
distribution has been obtained in |4]] and reads 

X (k)= J dbj dE J dP X (E, P; b) 6(yJ(E 2 - P 2 )/s - k) 6(VE 2 - P 2 - M min ) . (16) 



Then, by fitting the only existing x( K ) data [|15[ at ^fs = 16.5 GeV, we fix the parameter 



a of the model as a = 21.35. A comparison with the data is shown in Fig.l, where we 
have also put the result of 0. It is seen that the original version of IGM already gives a 
reasonable description of the data, but the inclusion of the impact-parameter fluctuation 
drastically improves the agreement. The latter enhances the small-K events and makes the 



3 For a recent comparative study of the main existing models on inelasticity, see for example 



overall shape flatter. The enhancement of large-K events is simply due to the larger value of 
a which is necessary for an overall fitting now. 

The inelasticity distribution x{k) has been obtained in || and reads 

X (k) = JdbjdEj dP X (E,P;b)5((E + P)/^7s-k)6(VEi-Pi-M min ). (17) 

Often, the inelasticity is defined in the lab. frame but, except when k — > 1, the difference 
between k defined in this frame and the one given in the e.v. frame is quite negligible. So, we 
will not make any distinction here and compute everything in the latter. We show, in Fig.2, 
the results for several pA collisions at y/s = 550 GeV. No accelerator data at such a high 
energy exists, but it is seen that x{k) is nearly k independent for pp, in agreement with ISR 
data [16| , 17] . In a recent cosmic-ray experiment [|18[], hadron-P6 inelasticity distribution 
at an average energy of (y/s) = 550 GeV has been estimated. The result^ is x(k) — 
(-0.25 ± 0.50)(1 - A;) L85±L77 + (3.1 ± 1.7)£; L85±L77 . We find that the qualitative features of 
our result agree with this estimate, except in the low-k region. Some of the origins of the 
discrepancy may be the difference between it — Pb and p — Pb collisions and the absence of 
the leading-particle fragmentation in our model. However, the uncertainty in their estimate 
is quite large. 

We show in Fig. 3 the average inelasticity (k) as function of -^/i, for several target nuclei. 
In Ref . [[| , several models for inelasticity have been compared with their estimates obtained 
with cosmic-ray data, which show slowly decreasing behavior with E . According to their 
comparison, IGM is the model which predicts the most quick fall of the average inelasticity 
with E Q in clear conflict with their estimates. Our curves in Fig. 3 still decrease as y/s 
increases but, compared with the results of ||, the energy-dependence is quite small now 
and compatible with the estimates obtained in ||. The main origin of this contrast is the 
factor a™ el which has been dropped out in (|7]), because it is not necessary in our version. 



4 In H|, we cited the preliminary result communicated by E. Shibuya, a member of the collabora- 
tion, which is slightly different. 



B. Leading-particle spectra 



A related quantity is the leading-particle spectrum, as shown in Fig. 4 at yfs = 14 GeV 
. Since data on dependence are scarce, we have assumed an approximate factorization 
of xi{= 2pi/y/s) and pr dependences, 

E^a/dp 3 )^ f{x l )h{p T ) 1 (18) 

where 

f(Pl) = JdbJdPjdE X (E, P; b) 6(VE* - P 2 - M min ) 5( [ft - (E + P)} /2 - Pl ) , (19) 
and parametrized h{p T ) as 

h(p T ) = ((3 2 /2n) e-P* , (20) 



determining the average {3 by using the p T dependence of the data |0| . The curves obtained 
with these (3 values (with an interpolation for Al and Ag) are shown in Fig. 4. The result of 
H for pp is also shown for comparison. Again, it is seen that a reasonable agreement with 
the data is obtained with IGM only, but the inclusion of the impact-parameter fluctuation 
improves it remarkably. We did not put their curves for the other targets, but the behavior is 
similar, namely they are more bent showing a definite deviation from the data in the largest- 
x\ region. This is a consequence of the neglect of the peripheral events there. Some authors 



TT| , p!3|JT4]| have obtained good fits to pA data, but in those works it is not clear which is the 
connection to other relevant quantities such as momentum distributions, correlations,... of 
the secondary particles. Also, pp is usually treated as a separate case. 



C. Rapidity distributions 

Let us now study the effects of fluctuations we have just introduced on the final particle 
spectra. As mentioned in the Introduction, the so-called inclusive rapidity (or pseudo- 
rapidity) distributions are the averages over such event-by-event fluctuations for a given set 



of experimental initial conditions. However, in usual computations of these distributions, 
say, by use of hydrodynamic models, instead of taking such averages 

l^j = J db J dP J dE^(E,P;b) X (E, P; b) 9{y/{E* - P*)/s - n min ) , (21) 

one considers some average initial conditions and calculates the distributions, starting from 
them. Namely, 

f ((E>,<P> = 0,<S». (22, 

It is evident that only under very special conditions that these two quantities can coincide. 

In [JD, we studied the deviation of the latter from the more realistic former distribution 
(dN/dy), by adopting the one- dimensional Landau's hydrodynamic model for an ideal gas 
in order to compute the rapidity distribution dN(E, P, b)/dy for a definite initial conditions. 
Despite all the simplifications, this model is known to reproduce the main features of the 
measured momentum (or rapidity) distributions and has advantage of having an analytical 



solution over the whole rapidity range ||20|| . The only inputs of the model are the total 
energy, momentum and the geometrical size of the initial fireballs. Remark that we are 
talking about the central fireballs and not about the rest of the system. In the case the 
incident particles are nucleons, the latter appears most frequently as leading particles. In 
order to avoid additional complexities, let us consider only this case, namely pp collisions. 
The invariant momentum distribution of produced particles in a hydrodynamic model is 



usually given by Cooper-Frye formula El 



dN r 

e W.™ /(pW<k ' (23) 

where cr(Td) is a constant-temperature freeze-out hypersurface, fip^u^) is Bose-Einstein (or 
Fermi-Dirac) distribution, p M is the 4-momentum of the emitted particle and is the 4- 



velocity of the fluid. Although it is possible to use more realistic freeze-out criteria p2]-|25 



here we limit ourselves to the simplest choice (|23|) without sophistication. This is enough 
for our present purpose of studying how the initial condition fluctuations affect the final 
particle spectra. 



cx(Trf) in is determined by solving the hydrodynamic equation 

d^ v = , with = {e + p)u»u u - pg^, (24) 

with an appropriate equation of state (p = e/3 in our case) once the dissociation temperature 
Td is given. As for the "initial volume" for a fireball of mass M, in the lack of a better 
justification founded on a physical basis, we adopt in the present work 

K„ = ^, (25) 



which has been suggested by a phenomenological analysis |26j of the M dependence of 
average multiplicity data |27[] and also consistent with the M dependence of the momentum 
distribution data The initial temperature T is then computed by putting M into this 
volume. As remarked in nowadays we know that neither the hypothesis of instantaneous 
thermalization nor the appearance of extremely high values of the initial temperature are 
physically reasonable. However, in spite of these rather non- conventional initial conditions, 
many of the qualitative and the quantitative results (average multiplicity, particle ratios, 
momentum distributions, • ■ ■) are surprisingly good when compared with data. In our point 
of view, perhaps the equilibrium is attained at a later time when the system has already 
suffered some expansion, but then the temperature and the rapidity distributions at the 
onset of the hydrodynamic regime would be approximately those of Landau's model whose 
initial conditions correspond to high temperature and energy density if extrapolated back in 
time. So, for any practical purpose, we can use Landau's solution to describe the system. We 
emphasize, however, that the fluctuation effects which are the central object of the present 
study do not depend sensibly on such a choice. 

We show in Fig. 5 a comparison of the results obtained in this way for < dN/drj > 
and dN/drj at y/s = 53 GeV. It is seen that the rapidity or equivalently pseudo-rapidity 
distributions are very sensitive to the fluctuations in the initial conditions. The peak, in the 
case of dN/dn computed with one fireball of mass < M >, corresponds to the simple- wave 
solution. When the fluctuations are taken into account, such a peak is completely smoothed 



away. They also cause a widening and a lowering of the distributions. Although the main 
purpose of this work is just to show the influences of the fluctuations, we may also compare 



them with some data ||29|| . We see that the behavior of the first one is more similar to the 
data than the other one and the presence of the simple-wave peaks in each event does not 
invalidate the overall agreement with data. 

IV. MULTIPLICITY FLUCTUATION 

We have shown in the last two sections how the impact-parameter and the energy- 
momentum fluctuations in the initial stage of the collision affect some of the observables. 
However, there are many more quantities whose description cannot be given only in terms 
of the fluctuations considered up to this point. Even after the mass of the fireball has been 
defined, quantities such as the multiplicity, particle species, their momentum distributions, 
... vary from event to event and hydrodynamic model we have used in the last section only 
describes the average behavior. Under certain conditions (constant dissociation temperature 
T d ), it does give the moments of the multiplicity distribution |3(|, so in principle also the 



multiplicity distribution itself [J3TJ] , but not the fluctuating events. 

Let us discuss in this section how the multiplicity fluctuation may be implemented. One 
way of doing this is to conveniently parametrizing the multiplicity distribution for a fixed 
mass M and determining the parameters by imposing certain constraints. We choose a very 
simple parametrization for the multiplicity distribution 

tfj{M,z) = Az u e~ az , (26) 

where, as usual, 

i/>(M,z) = (N(M)) P N (M), (27) 
z — N/ (N(M)) 

with Pn(M) indicating the probability of a fireball of mass M decaying into N charged 
particles, and impose the conditions 



1 f°° 

Y,Pn{M) = 1 --/ ^(z)dz = l, (28) 



N 2 Jo 



J2NP N (M) = (N(M)) --/ s^O*)^ =<*> = !, 



By substituting fl26|) into these equations, we obtain 

i/ = a-l, (30) 

A =^o- (31) 

The moments (z n ) = (N(M) n ) / (N(M)) n have been calculated in [31] and, in particular, 
it is found that (z 2 ) = 1 + a 2 / (JV(Af)), with a 2 = 1-105. What we can do is, once M is 
fixed, to produce events following this distribution by the use of the Monte Carlo method. 
In doing so, we have indeed to consider also the charge fluctuation. 

The overall multiplicity distributions calculated in this way is shown in Fig. 6. As seen, the 
results reproduce quite well the qualitative features of the data in all the ISR energy region. 
They are slightly narrower than the data and, as the energy increases, the discrepancy 
becomes more pronounced, but probably the data begins to suffer the influence of the mini- 
jets there. 



V. ENERGY AND MOMENTUM FLUCTUATION OF THE OBSERVED 
PARTICLES: PHASE-SPACE MODEL 

The multiplicity fluctuation, discussed in the preceding section, does not manifest it- 
self only in the multiplicity distribution. With the inclusion of this fluctuation, we are 
considering that the momentum distribution of the secondary particles for a given M is a 
superposition of distributions with different multiplicity N. Moreover, even with a fixed N, 
the momentum distribution will vary from event to event. As will be shown below, there 
are observables such as the forward-backward multiplicity correlation, which depends on 



this kind of fluctuation. But, then we faces the following problem: "How to compute the 
rapidity distribution of a system having a definite mass M and a definite multiplicity iV 
and in an event- dependent way?" This question did not arise when computing the inclusive 
distribution, because the hydrodynamic model does take such a fluctuation into account, 
as mentioned in the preceding section. What it does not do is to generate each fluctuating 
event. 

We propose to use the one- dimensional phase-space model to generate these events. 
First, one- dimensional because we know from the data that the momentum distributions 



in high-energy is essentially longitudinal. In a previous work [33|, we have shown that the 
rapidity distributions predicted by the one-dimensional phase-space model, given M and 
N, are approximately Gaussian, as in hydrodynamic model. It was also shown that these 
are insensitive to a certain class of dynamical factors introduced in the model. Although 
it is indeed not guaranteed a priori that the superposition of these distributions shall give 
the one obtained by the hydrodynamic model, the result of this model has the qualitative 
features of that one and an advantage of being a sum of event-dependent distributions with 
fixed M and N. 

Given a mass M and a multiplicity N, the one-dimensional phase-space model tells us 
that the probability of finding an event with the particles in the longitudinal momentum 
intervals 

\pi,Pi + dpi), i = l,---N, (32) 

is given by 

" Np -i^kk---^% p ' )s i E '- M) ' (33) 

where 

The invariant one-particle distribution is then given by 



dN = dN = NRn.^M') 

dp dy R N (M) ' 1 j 



normalized to iV and where M' = y(M — E) 2 — p 2 represents the invariant mass of the 
system after subtracting the observed particle. The correct expression of the probability 
would be ( |3"3"D with some dynamical factor. In [03|, in order to simulate the hydrodynamic 



motion, we have included a factor f(y) = a.e~^ mT<x ' slLy for each particle, where y is the 
rapidity, and shown that this is entirely irrelevant, the final result being the same. 

In Fig. 7, we show one-particle pseudo-rapidity distributions in pp collision at 53 GeV 
given by (0). Since we have completely neglected any dynamical factor, we could not 
expect to obtain a perfect agreement with the data. However, it is seen that the qualitative 
features of the data are reproduced. 

A. Forward-backward correlation 

One of the data, which cannot be understood without the multiplicity fluctuation dis- 
cussed in Section IV and nicely reproduced with the longitudinal phase-space model proposed 
here, refers to the so called forward-backward charged multiplicity correlation. First of all, 
the data show that the charged multiplicities in the two hemispheres are very little corre- 
lated jTjJ , presenting a large fluctuation. The correlation is usually presented as a graph of 



the average charged multiplicity in the backward hemisphere < iV& > as a function of the 
effectively observed charged multiplicity Nf in the forward hemisphere as in Fig.8. 

To begin with, remark that in the usual application of the hydrodynamic model, without 
any fluctuation taken into account, the graph would reduce to a single point. The consider- 
ation of the impact-parameter and the initial-state fluctuations, implemented through IGM 
as we are proposing, improves considerably the agreement with the data. Notice that the 
momentum fluctuation of the fireball is essential, otherwise the correlation would be too 
large. Now, without the inclusion of the multiplicity fluctuation and especially the momen- 
tum fluctuation of the final particles, the correlation begins to deviate from the experimental 



trends for large values of Nf , because, in that case, large Nf means large M with the fireball 
sitting more or less in the center of mass, with a symmetrical distribution of particles. 



VI. CONCLUSIONS AND FURTHER OUTLOOKS 

We have investigated, in this paper, effects of several kinds of fluctuations which appear 
in hadronic and nuclear collision, analyzing different kinds of observables. The Interacting 
Gluon Model, improved by the inclusion of the impact-parameter fluctuation and comple- 
mented by an appropriate hydrodynamic model seems to describe well the the bulk of the 
phenomenon, such as the cross-section, inelasticity, leading-particle spectrum, average mul- 
tiplicity and the inclusive momentum distribution. Other quantities depend explicitly on the 
multiplicity and the final-particle momentum fluctuations. In this paper, we have treated 
the first one by considering the thermodynamics of the fireball and then parametrizing the 
multiplicity distribution for each mass M in a convenient way. More microscopic descrip- 
tion of this fluctuation would also be possible. We have treated the momentum fluctuation 
of the final particles by using the one-dimensional phase-space model. Although not com- 
pletely satisfactory, this tretement has shown the importance of such a fluctuation in an 
event-dependent basis. 

There are many other properties which clearly depend on fluctuations. One of these 
quantities is the so-called semi-inclusive rapidity distribution, namely, distribution with a 
fixed multiplicity interval. In principle, it would be possible to obtain this distribution with 
the ingredients we have considered here, but there is something which is missing. Especially 
in the low multiplicity intervals, the effects of the diffractive processes cannot be neglected. 
We are studying how to incorporate the diffractive processes in IGM, in a consistent way. 
Also, to compare with the data, a treatment which is somewhat more realistic than ( p^f ) is 
required. Another interesting quantity that could be studied, which certainly depends on 
event-by-event fluctuations we have discussed, is the Bose-Einstein correlations of produced 
particles, or the so-called Hanbury-Brown Twiss effect f35|, frequently used in heavy- ion 



collisions to infer about the space-time developpment of the hadronic matter which is formed 
in such collisions. Suggestions for such a study has been given by M. Namiki et al. in 0. 

In this paper, we have completely neglected the fragmentation of the leading particles, 
which certainly give non-negligible contributions to the semi-inclusive rapidity distributions 
mentioned above, especially when the multiplicity is small. In the case of nucleus-nucleus 
collisions, how to treat the fragmentation of the leading nuclei is a completely open question. 
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Figure Captions 



Fig.l: ^-distribution for pp at y/s = 16.5 GeV. The data are from |15 l. The solid line is 
our result, whereas the dashed one is from [|]. 

Fig. 2: Inelasticity distribution for pA collisions with several targets at -^i = 550 GeV. 

Fig. 3: Energy dependence of the average inelasticity for pA collisions. 



Fig. 4: Leading-particle spectra as fuction of xi at pt = .3 GeV. The data are from [19 
at \/s = 14 GeV. The solid curves are our results, whereas the dashed one is from 
The slope parameter has been extracted from |19[ as (3 = 4.20, 3.22, 3.21, 3.26, 3.47 
and 3.78GeV _1 for p, C, Al, Cu, Ag and Pb targets, respectively. 

Fig. 5: Pseudo-rapidity distributions calculated in the usual procedure (solid line) and 
with fluctuations (dashed line) at = 53 GeV. Experimental data |29j are shown for 
comparison. 

Fig. 6: Charged-particle multiplicity distribution in pp collisions at yfs = 30 and 62 GeV, 
compared with data |32[ . 



Fig. 7: One-particle inclusive pseudo-rapidity distributions for pp collisions at y/s = 53 
GeV, computed by using the longitudinal phase-space model. The solid curve repre- 
sents the distribution without any fluctuation, whereas the dashed one is the result 
with all the fluctuations included in the way described in the text. The data are from 
PH. 



Fig. 8: Forward-backward multiplicity correlation in pp collisions at yfs = 24 GeV, com- 
puted with the longitudinal phase-space model (□), compared with the result without 
the multiplicity and final-particle-momentum fluctuations (*). The data (•) are from 
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